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Abstract. Equipping the tangent bundle TQ of a manifold with a symplectic 
form coming from a regular Lagrangian L, we explore how to obtain a Poisson- 
Nijenhuis structure from a given type (1, 1) tensor field J on Q. It is argued that 
the complete lift J'^ of J is not the natural candidate for a Nijenhuis tensor on TQ, 
but plays a crucial role in the construction of a different tensor i2, which appears to 
be the pullback under the Legendre transform of the lift of J to T*Q. We show how 
this tangent bundle view brings new insights and is capable also of producing all 
important results which are known from previous studies on the cotangent bundle, 
in the case that Q is equipped with a Riemannian metric. The present approach 
further paves the way for future generalizations. 



1 Introduction 

There is a well established theory of bi-Hamiltonian systems on the cotangent bundle 
T*Q of a manifold Q or, more generally, on a symplectic or Poisson manifold. In par- 
ticular, there is a link between separability of the Hamilton- Jacobi equation and cer- 
tain classes of bi-Hamiltonian or quasi-bi-Hamiltonian systems, in which more specif- 
ically Poisson-Nijenhuis structures ^H] play a prominent role. The immediate source 
of inspiration for the present work is a series of recent papers in this general field in 
which relations have been explored between such things as: bi-differential calculi, com- 
plete integrability, Stackel systems, compatible Poisson structures on an extended space, 
Gelfand-Zakharevich systems, so-called special conformal Killing tensors, and so on. See 
[3llinilHllllinilIIllI2llIllliniinil21l2nifora non-exhaustive list of recent contributions. 
In the more direct applications of these theoretical developments, the dominant geomet- 
rical space is a cotangent bundle, which of course comes equipped with its canonical 
symplectic or Poisson structure, and on which a compatible Poisson structure is obtained 
via the lift J of a type (1, 1) tensor field J on the base manifold, thus giving birth to a 
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Poisson-Nijenhuis structure. The link with a class of Stackel systems requires the avail- 
ability of a metric g on Q, with respect to which J has the property of being a special 
conformal Killing tensor (a concept which is intimately linked with what is called Benenti 
tensor after the extensive work of Benenti on Hamilton- Jacobi separability (see e.g. 

mm))- 

The point we want to emphasize now, however, is that some of these applications clearly 
come from meaningful questions about dynamical systems living on a tangent bundle TQ. 
This is for example the case with the theory of cofactor pair systems, as developed origi- 
nally by Lundmark on a Euclidean space [T71[TH] and generalised to (pseudo-)Riemannian 
manifolds and in more geometrical terms by Crampin and Sarlet jHj. The physical back- 
ground for these systems is a kinetic energy type Lagrangian on TQ for which (at least 
in one possible interpretation) admissible non-conservative forces are being sought, in the 
sense that the resulting Newtonian system admits two quadratic first integrals, which in 
turn can generate a whole family of integrals in involution. The cofactors of the Killing 
tensors coming from these integrals then determine the special conformal Killing tensors 
which give rise to Poisson-Nijenhuis structures; we are then looking at examples of so 
called bi-quasi-Hamiltonian systems [H]. 

A type (1, 1) tensor field J on Q also has a natural lift to TQ] it is usually called the 
complete lift and we will denote it by J^. The difference with T*Q of course is that TQ 
does not carry a canonical Poisson structure. However, a symplectic form is available (and 
can be constructed by pure tangent bundle techniques) as soon as a regular Lagrangian 
is given, which could for example be the kinetic energy Lagrangian coming from a metric 
on Q. It does, therefore, perfectly make sense to explore the possibility of obtaining 
Poisson-Nijenhuis structures on TQ by natural tangent bundle constructions. One of the 
goals we have in mind for the future is to arrive at a generalization of the theory of special 
conformal Killing tensors from Riemannian to Finsler spaces. The primary objective of 
this paper, however, is to set the ground for future developments by trying to understand 
in detail how the results one is by now familiar with in a cotangent bundle environment, 
can be obtained in a natural way by pure tangent bundle techniques. We shall see that 
this different way of approaching the subject offers new insights anyway. In fact, some of 
the preliminary considerations lead to results which are valid for arbitrary Lagrangians, 
not just kinetic-energy type ones. In this respect, we are to some extent joining the 
interest in Poisson structures on a tangent bundle which is present also in recent work by 
Vaisman [2111301 Ell- 
in Section 2, starting from a given Lagrangian L on TQ and a type (1, 1) tensor field 
J on Q, we show that J^S determines an alternative almost tangent structure and use 
this to construct another type (1, 1) tensor R on TQ. The full characterization of R 
is developed in Section 3, making use of the connection provided by the second-order 
equation field P, coming from the Lagrange equations. The eigenspace structure of R 
is discussed in Appendix B, whereas further properties of general interest are derived in 
Section 4. In Section 5, we specialize to the particular case of a Riemannian manifold 
and the associated kinetic energy Lagrangian and explore how various known properties 
make their appearance within such a tangent bundle approach. Some unexpected new 
features come forward which are further discussed in Appendix A. The road map to future 
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developments is sketched in the final section. 



2 Generalities 

Suppose we are given a regular Lagrangian L on TQ and a type (1, 1) tensor field J 
on Q. L comes with its associated symplectic form uol, and J determines a tensor field 
J'^, its complete lift, on TQ. One may wonder whether these data can give rise, under 
some circumstances, to a compatible Poisson structure. For that, J'^ should be symmetric 
with respect to ul, the so-called Magri-Morosi concomitant must vanish (see |23 IT9]). 
and also the Nijenhuis torsion of must be zero (which is equivalent to the torsion of 
J being zero). It seems to us, however, that this is not the most interesting path to 
pursue. Experience in a variety of apphcations (see for example [71 I2H1) has shown 
that interesting type (1, 1) tensor fields i? on a symplectic manifold (M, a;) arise from the 
construction of a second 2-form Ui and the determining formula: 

= i^ui, G XiTQ). (1) 

A direct advantage of such tensor fields is for example that they automatically have the 
required symmetry property 

ujiR^,ri)=uj{^,Rr]). (2) 

Furthermore, vanishing of the Magri-Morosi concomitant then is equivalent to dui = 
|inj . after which we are left with the condition = 0. Although the canonical lift J 
of J to T*Q can be defined (see |Zj) via a relation like ((H), this does not seem to be the 
case for J'^ on TQ. We therefore start our investigation with an exploration of possible 
natural ways of constructing a second closed 2-form ui from the given data on TQ. 

For completeness, we list a number of useful properties of J'^, starting with defining 
relations with respect to the action on complete and vertical lifts of vector fields on Q: 
for X,Y eXiQ), 

r{x') = {jxy, rix"") = {jxy. (3) 

Using the bracket relations 

[X\ = 0, [X^ Y^] = [X, r]^ [X^ Y'^] = [X, YY, 
it easily follows that 

Cx^r = {CxJy, Cxvr = {CxJf. (4) 

Also, for the Nijenhuis torsion we have, 

iVj.(X^,F^) = 0, 
Ar^c(X^n = {Nj{X,Y)r, 
NAX\Y^) = NAX^Y^) = {Nj{X,Y))\ 

from which it follows that Njc = <^=^ Nj = 0. 
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It is imperative to relate J*^ to the canonical type (1,1) tensor field S on TQ (the so-called 
vertical endomorphism) , which satisfies S{X'^) = , S{X^) = 0. It is easy to see that 
J'^ commutes with 5", in the sense of endomorphisms on X{TQ), but also in the sense of 
the Nijenhuis bracket: 

[ J^ s] = 0. (5) 

It follows that also the corresponding degree 1 derivations commute, meaning that: 

dsdjc = -djc ds- (6) 

Finally, it is easy to verify that for any J, 

Njcs = 0. (7) 

In fact, we can make the following more complete statement in that respect, which is 
trivial to prove, and essentially says that J'^S very much has the same properties as S. 

Lemma 1. We have (J'^S)'^ = and Njcg = 0. Furthermore, if J is non-singular, J'^S 
determines an integrable almost tangent structure. 



Now recall the role which the vertical endomorphism plays in the definition of the Poincare- 
Cartan 2-form cul' we have, 

ujl = d{S{dL)) = ddsL. 

[We make no notational distinction between the action of a type (1,1) tensor field on 
vector fields and its dual action on 1-forms, but one should keep in mind that the order 
of composition of such action changes in passing from one interpretation to the other.] It 
then looks perfectly natural, given J and the sort of alternative integrable almost tangent 
structure which it creates, to consider the closed 2- from uji, defined (with various ways of 
writing the same expression) by, 

uji = d{Sr{dL)) = d{S{djcL)) = d{r{dsL)) = ddjcgL. (8) 

And so, with L and J as data, the type (1, 1) tensor field R which will carry our attention 
is defined by 

iR{^)ddsL = i^ddjcgL. (9) 

We know that it will define a Poisson- Nijenhuis structure if and only if — 0. 

The first objective now must be to obtain a reasonably practical description of i?, for 
example by recognising its action on complete and vertical lifts. In fact, we believe that 
it is better for general purposes, to make use of horizontal and vertical lifts, rather than 
complete and vertical lifts. For that, of course, one needs a connection, but there is one 
available, namely the non-linear connection associated to the Euler- Lagrange equations 
of L (being second-order differential equations on TQ) . 
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3 Making use of a connection 



As is well known, every second-order equation field 

in particular the one coming from the regular Lagrangian L, defines a horizontal distri- 
bution, with connection coefficients 



dp 

p. = ill) 



As a result, every vector field ^ on TQ has a unique decomposition of the form ^ = 
X" + y^, where X, Y are vector fields along the tangent bundle projection r : TQ Q. 
The C°°(TQ) module of such fields is denoted by X{t). An extensive calculus along the 
projection r was developed in We recall here some basic features of this calculus, 

which will be needed in what follows. 

Interesting derivations and tensorial objects along r are discovered by looking at the 
decomposition of Lie brackets of vector fields on TQ. We have, for example, that 

[X\Y^] = (12) 
[X^,y^] = {D^Yf-{D^X)", (13) 

= {[x,Y]^r+n{x,Yr. (m) 

Here, TZ is the curvature of the non-linear connection, a vector valued 2-form along r, 
and are the horizontal and vertical covariant derivative operators, which act on 
functions F G C°°{TQ) as 

D5F = X^(F), D^F = X"'(F). (15) 

The horizontal and vertical brackets of elements of X{t) (in the present situation of a 
connection which has no torsion) are given by 

[X, y], = D^F - D^X, [X, Y]^ = Djr - D?X. (16) 

Other brackets of interest are 

[r, x^] = -X" + (vx)^, [r, x^] = (vx)" + $(x)''. (17) 

Here $, a type (1, 1) tensor along r, is called the Jacobi endomorphism, and V is the 
dynamical covariant derivative, which on functions acts like F. 

There exist lift operations on many other tensorial objects. We mention two more con- 
structions of interest now. First there are the horizontal and vertical lifts of a type (1,1) 
tensor field U along r, determined by 

[/^(X^) = [/(X)^, U"{X") = U{XY, (18) 
U^{X^-) = 0, f/^(X^) = U{XY. (19) 
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Next, if is a symmetric type (0, 2) tensor field along r, its Kahler lift is a 2-form on TQ 
determined by 

g-{X",Y")=g-{X\Y^) = Q, (20) 
gK^X\ Y") = -g'^iX", Y^ = giX, Y). (21) 

In fact, the Poincare-Cartan form ujl is precisely the Kahler lift of the Hessian of L, 
defined intrinsically by 

g = D^'D^L, (22) 

where is the vertical covariant differential defined (on any tensor T) by D^T(X, . . .) = 
Dj^T(. . .). Observe that the complete lift J'^ of a (1, 1) tensor J on Q, can be written as 

J- = J^ + {VJf. (23) 

Other, more specific properties of interest in the calculus along r will be recalled when 
appropriate, but we should at least refer here also to the existence of a canonical vector 
field along r, the total time derivative operator T = u^d/dq^, whose vertical lift is the 
Liouville vector field A, whereas its horizontal lift is a second-order equation field, which 
need not be the one we started from; the two coincide for a quadratic spray. 

We can now start the computation of the structure of the tensor field R defined by (jHI). 
To begin with, observe that from (fTR|) . (fT!Hl and (jSBI), it easily follows that 

J^(X^) = {JXy, (24) 
r{X") = {JXY + VJ{XY. (25) 

Using the standard notation 9l for c/^L, we have 

0^(X^) = 0, eL{X") = dL{X^) = X^{L) = D^L = d^L(X). (26) 

Note that in the second of these relations, a computation on TQ in the end is replaced by 
one involving a vector field and 1-form along r. The point is that 6*/,, being semi-basic, 
can be regarded as a 1-form along r as well; the defining relation then reads 9l = d^L, 
where d^ is the vertical exterior derivative. The latter is completely determined by the 
following action on functions F and 1-forms a along r: 

d^F(X) = D^F, rf^a(X, Y) = D^a(r) - D^a{X). (27) 

Similar relations hold for the horizontal exterior derivative d" . The action of these exterior 
derivatives extends to vector- valued forms as well; it suffices for our purposes to know that 
for d^ or d" acting on a type (1, 1) tensor field along r, the defining relation is formally 
the same as in the second of equations 

It is easy to see that ui = d^J'^OL) gives zero when evaluated on two vertical vector fields. 
Next we have, passing as before from 6*^, regarded as 1-form on TQ, to its interpretation 
as 1-form along r, 

u,{X\Y-) = Cxv(dr.{{JYY))-eL[j\[X\Y"])) 

= D^(^^(jr)) - eL(j(D-r)) = D-e^(jr), 
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from which it follows in view of ()22|) that 

uj^{X\Y")=g{X,JY). {2i 
Proceeding in the same way, we get 

= B"^{9,iJY)) -B-(9LiJX)) -9l[j{[X,YI 

= D%iJ9L{Y)) - D?(J^^i(X)) - ^^l(j(D5F - D?X)^ 
= B"AJ0L)iY)~D-iJ9,)iX). 

It follows that 

u,iX-,Y") = d"iJeL)iX,Y). (29) 
Proposition 1. The type (1, 1) tensor field R, defined by ^ has the following structure: 

/?(X^) = (JX)'' (30) 
R{X") = {JX)" + {UXy, (31) 

where J is the transpose of J with respect to g = D^D^L, i.e. g{JX, Y) = g{X, JY), and 
U is the tensor field along r, determined by 

g{UX,Y)=d"{jeL){X,Y). (32) 

Proof. It is sufficient to take horizontal and vertical lifts of basic vector fields for finding 
the tensorial structure of R. We have ujl{R{X'^),Y^) = and 

coLiRiXn,Yn=co,iX^,Yn=giJY,X) = g{Y,JX)=g''{(jX)\Y"), 

from which ()30|1 follows. Likewise 

ujL{R{Xn,Yn = -g{JX,Y)=g-{{JXr,Yn, 

from which it follows that R{X") = (JX)" + {UXY , for some U. Subsequently, using 
iPH) and IPU)). 

UL{R{X"),Y^)=UL{{UX)\Y")=g{UX,Y) = d"{jeL){X,Y), 

which completes the proof. □ 

Note that it follows from the skew-symmetry of the right-hand side in (I32|l . that U = —U. 

In Appendix B, we investigate the eigenspace structure of R and the explicit construction 
of so-called Darboux-Nijenhuis coordinates. 
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4 Further properties of the tensor field R 



Proposition 2. We have R = J'^ <^=^ J = J and U = V J. 

Proof. The result follows immediately from comparison of ()30II31|) with ()24II25|) . □ 

A natural question which arises is whether R, in general, could commute with S, just as 
J'^ does, either in the algebraic sense or with respect to the Nijenhuis bracket. 

Proposition 3. RS = SR J = 1. 

Proof. Using S{X^) = and S{X") = X^, the result follows immediately from the 
characterization of R in Proposition 1. □ 

Recall that the Nijenhuis bracket is defined by 

[R,S]{^,r]) = [R^,Sv] + [S^,Rv] + {RS + SRm,v]) 

- Rm, f]] + [e, Sf]]) - V] + [e, Rv])- 

Proposition 4. Assuming J = J, so that RS = SR, we have [R, 5*] = <^=^ = d"J . 

Proof. That [i?, S] vanishes on two vertical lifts is trivial. Again, it suffices for such 
calculations to consider lifts of basic vector fields (vector fields on Q) , rather than vector 
fields along r. Since J is basic as well, one then easily verifies, making use of the bracket 
relations ^ and (O, that also [R, S]{X'' ,Y") = 0. Next, we have 

[R,s]{x",Y") = [(jx)^ + (f/x)'',r^] + [x^,(jr)^ + (t/r)''] 

+ 2RS([X, F]^^) - Y"] + [X", y^]) 

- s([ijx)" + {uxy, Y^] + [X", {JYY + [UYY]) . 

Using the bracket relations ()12lll4p . plus the fact that J, and by assumption also X en 
F, are basic, this readily reduces to 

[R,S]{XYY") = (D-t/(r) -D-^(X))" + (J(D5r-D-X))^ + (D-(JX) -D5(Jr))^ 
which in view of properties such as ^I7\ for type (1,1) tensors, can be written as 

[R, S]{X", Y") = {dm{X, Y) - d"J{X, Y)Y. 

The result now follows. □ 

Another obvious question is under which circumstances R is truly a recursion operator for 
symmetries of F. For that, we compute C^R- Taking the Lie derivative with respect to 
F of the defining relations (j3(JH3ip and making use also of the properties (fTTj) . it is fairly 
straightforward to verify that 

£ri?(X^) = {J-J)iXf + iU + VJ){X)\ (33) 
CrR{X") = {VJ -U){X)" + {VU + (^J -M){XY. (34) 

The following, therefore, is an interesting immediate result. 
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Proposition 5. CrR = J = J, f/ = VJ = 0, <I>J = J$. □ 

Next, we address the question of recognizing the conditions under which R has zero 
Nijenhuis torsion. A direct computation of Nji, through its action on horizontal and 
vertical hfts is extremely tedious and therefore not worth the effort, since we can actually 
rely on existing results concerning the complete lift J of J to the cotangent bundle T*Q. 
Let Leg : TQ T*Q denote the Legendre transform coming from the regular Lagrangian 
L. 

Proposition 6. The tensor field R defined by (0j is directly related to the complete lift 
J on T*Q: as a matter of fact we have R = Leg* J. 

Proof. The defining relation of J, as described in reads (where ^ this time denotes 
an arbitrary vector field on T*Q), 

ij^dO = i^CjvdO, (35) 

and clearly has the same structure as 0- Here, is the vertical lift of J, which is a 
vector field on T*Q (cf. ^^^1 6 of course is the canonical 1-form on T*Q. It is easy to 
see from the coordinate expression that ijvd6 = J6, so that (jH^ implies ij^dO = i^dJ6. 

But it is equally trivial to verify in coordinates that Leg* J 9 = J'^O^ = djcgL. The result 
then immediately follows from taking the puUback under Leg of the new representation 
of dSl. □ 

As an immediate consequence of the fact that Nj = -v^ Nj = 0, which incidentally 
requires also a fairly tedious calculation (cf. [7 ), we now come to the following conclusion. 

Proposition 7. A^^ = Nj = 0. □ 

It is of some interest to recall here the following characterization of Nj = 0. 
Lemma 2. Nj = if and only if for all basic vector fields X,Y, we have 

^jxJiY) - JmJ{Y)) = B"jyJ{X) - J(D-J(X)). (36) 

Proof. Taking into account that [X, Y] = [X, Y]^ for basic vector fields, we have 

[JX, JY] + J2[x, Y] - J{[JX, Y] + [X, JY]) = B"jx{JY) - B"jy{JX) 
+ J\D-Y - D-X) - J(D-^r - D-(JX) + D- (jr) - DJ^(X)) 
= DJ^ J(F) - J(D- J(r)) - B"jyJ{X) + J(D- J(X)), 

which gives the desired result. □ 

To conclude this section, we look at an interesting case of non-vanishing C^R- Observe 
first that F, by construction, is the Hamiltonian vector field associated to the energy 
function = A(L) — L, with respect to the symplectic form ul = dO^, i.e. we have 

i^dOL = -dEi. (37) 
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A general theorem proved in JU] then imphes that in the present context, the fact that 
the 2-form uji is closed, is equivalent to stating that 

^CrRdOL = -'2ddREL. (38) 

Since R is symmetric with respect to ul, the same is true for C^R, so that invariance of 
R is equivalent to having ddnE]^ = 0. The latter was the starting point for an application 
of a bi-differential calculus in ^Oj, to which we shall return in the next section. The more 
general situation of a gauged bi-differential calculus in |10j, corresponds to the assumption 
that for some basic function /, 

ddnEL = df A dEi. (39) 

Via the equality this assumption is equivalent to stating that (cf. Prop. 5.3 in [lOj) 

CrR = r ® rf/ - e/ ® dEL, (40) 

where C,f is the Hamiltonian vector field associated to /. It easily follows, from is^^^ujl = 
—S{i^^.ujL) = S{df) = 0, that ^/ is vertical, say ^/ = Xf^ for some Xf along r. Then, 

z^M"^") = (7-(X/, F-) = g{Xj, Y) = -Y"{f). 

Hence, in terms of fields along the projection r, Xf is defined by 

g{Xf,Y) = -Y"U) = -d"f{Y) or Xf}g = ~d"f. (41) 

In the next section we specialize to the case that the Lagrangian comes from a (pseudo) 
Riemannian metric on Q, and will focus most of the attention on the characterization of 
so-called special conformal Killing tensors in their tangent bundle manifestation. 



5 The Riemannian case 

Let 5^ be a symmetric, non-singular type (0,2) tensor field on Q and put L = ^gijU^u^ . 
The non-linear connection defined by the Euler-Lagrange equations then is the (linear) 
Levi-Civita connection of g, i.e. the connection coefficients (fTT|) are of the form 

r} = ^},«^ (42) 

where F*^ are the classical Christoffel symbols. 

It is important to understand first how the fundamental covariant derivative operators 
of the calculus along r, as referred to in the previous section, relate to classical tensor 
calculus in this case. For example, if J is a type (1, 1) tensor field on Q, then both the 
dynamical covariant derivative VJ and horizontal covariant derivatives like DJJ relate 
to the classical covariant derivative Jj|^, as follows: 

(VJ)} = JJ|,«^ (D- J)j = j;^,X\ (43) 
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In fact, in this situation we have V = Df^. It follows that, in particular, Vg = and 
Dx9 = O5 {and of course also Dxd = because g is basic). Another specific property 
of the case of a quadratic spray is that the so-called deviation, which in the language of 
the calculus along r is VT, is zero. It is actually of interest to list all covariant derivatives 
of T here: 

VT = 0, D^T = X, DJT = 0. (44) 

These properties are easy to verify in coordinates, but let's take the opportunity to men- 
tion also the general commutator identity 

[V,D^]=D-^-D-, (45) 

which can be used to show in a coordinate free way that the first two relations (j44|l imply 
the third. 

We now look at the characterization of R in this context, more particularly the specifica- 
tion of the tensor field U. Note first that everything should now be easily expressible in 
terms of the metric since we have, for example, that 

eL{X) = g{T,X) or = Tj^?. (46) 

We should keep in mind also that L = El is a first integral, so that VL = T{L) = 0, and 
that it further follows from using the commutator property [V, ci^] = —d", that 

= V^L = Vd^'L = -d"L. (47) 

Now, concerning the determination of f/, we have 

d''{jeL){X,Y) = D5(J0z.)(n-D?(^^L)(X) 

= D5(^7(T, JY)) - g{T, J(D^r)) - D?(^7(T, JX)) + g{T, J(D^X)). 

Taking into account that T)^g = and D^T = 0, we conclude that U is determined by 

giUX, Y) = giT, J(F) - J(X)) = ^(T, d"JiX, Y)). (48) 

In coordinates, we of course work with the adapted frame of horizontal and vertical vector 
fields on TQ and their dual 1-forms, which are 

{^^ = |i - ^, = ^} ' {^^^ = + ri^^n • (49) 

The coordinate expression for R then becomes 

R = X.Hi ® dq^ + 7jV, ® r]^ + U;V, ® dq\ (50) 

where 

U; = 9'\Jk\j - JT\k)9miu'. (51) 

It is interesting to have intrinsic expressions also, which implicitly determine the vertical 
and horizontal covariant derivatives of U . 
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Proposition 8. For Z G X{t), D|f/ and Df f/ are determined by 

(?(D^[/(X),F) = g{Z,D"^JiY)-D^JiX)), (52) 
g{D"zU{X),Y) = g{T,D''D"J{Z,X,Y) -D"D" JiZ,Y,X)), (53) 

Proof. The proof is a straightforward computation, which starts from PHj) and takes 
the properties (jl^ into account, remembering further that 

D"D"J{Z, X, Y) = Df J(r) - B^H^J{Y). □ 

Let us come back in this case of particular interest to the invariance of R, or the more 
general assumption ()40j). If h = ^g^-'PiPj is the corresponding Hamiltonian, we know 
from the results of the preceding section that CrR = is equivalent to ddjh = 0. It 
was shown by a coordinate calculation in |T^ that the latter condition is equivalent to 
J = J and VJ = 0. This is somewhat surprising now, since Proposition 5, in general, 
imposes more conditions for having C^R = 0. Of course, if V J = 0, it follows from (j43j) 
that also DJJ = and thus from (pSj) that also U = 0. This implies in particular, from 
Proposition 2, that when R is invariant, it must be equal to J'^. But the more interesting 
information which follows from comparison with Proposition 5 is that apparently, when J 
is symmetric and parallel in the Riemannian case, it will automatically commute with the 
Jacobi endomorphism $. This is not a trivial property to recognize. We therefore propose 
to verify in Appendix A by a direct calculation that it is indeed a correct statement. That 
direct proof is of interest in its own right, because it illustrates how one can proceed with 
an integrability analysis in this context. 

Now assume that (jinjl holds for some / G C°°{Q) and where ^/ = Xf^ , with Xf defined 
by ()4H1 . In the present situation, we further have El = L and F = (since VT = 0). 

Theorem 1. Under the present circumstances, the tensor field CrR is of the form ^C\ ) 
if and only if J = J and further satisfies 

\/J =\{T®d"f ~Xf®dL). (54) 

In addition, U then is of the form 

U = -\{T®d"f + Xf®dL), (55) 

and J further has the property 

- J^ = \{T®Vd"f + VXj®eL). (56) 

Finally, the tensor field R itself then is given by 

R = r - A®df. (57) 

Proof. The right-hand side of (jlOJ, when evaluated on some X^ , results in —X^{EL)Xf^ = 
—9L{X)Xf^. Comparison with (jH^ shows that this requires J to be symmetric, plus the 
condition that 

U = -VJ -Xf® Ol. 
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Proceeding in the same way for an arbitrary horizontal argument X" , comparison with 
(jnH) reveals the requirements 

[/ = V J - T ® d"f, 

and 

Vt/ + $ J - J$ = -Xf (g) d"L = 0. 

Compatibility of the two expressions for U above, immediately leads to the conclusions 
and (jnSl)- The first of these puts a restriction on J, while the second in fact is then 
automatically satisfied. To see this, we compute (F\/J from Remember (see pij ) 

that for a vector valued 1-form along r, of the form a® X, an exterior derivative such as 
(F is computed via the rule: (F {a® X) = d^a® X - a ^ (FX. Now (Fd"f = -d"d'^f = 0, 
= / (the identity tensor), d^d^L = and finally also d'^Xj = from (gH). It follows 

that 

d"J = d'^VJ = -\d"f A /. (58) 

The defining relation ()48p for U then easily leads to ()55|) and can in fact also be rewritten 

as 

g{UX,Y) = ~\d''f ^OL{X,Y). (59) 

The final requirement that Vf/ should be equal to J$ — $J leads immediately, from 
to (jSSl), or can equivalently, from (jSHl), be expressed as 

($J - J$)J g = lVd"f A Ol. (60) 

The point is, however, that this again is not an extra condition, but a consequence of the 
fundamental condition (|^. To see this, recall that pUj) is equivalent to which in 
turn, when translated into the corresponding cotangent bundle property reads ddjh = 
df A dh. It was shown by a direct coordinate calculation in ^| that this condition 
requires the symmetric J to be a so-called special conformal Killing tensor. We shall 
verify in coordinates below that this is exactly the condition (j54p . Hence, (|56p must 
be a corollary and one could obtain it in a direct way by following the pattern of the 
integrability analysis in Appendix A. 

The final statement ()57p about R follows directly from comparison between ()24ll25p and 
knowing that J = J and using ^ and with = A. □ 

In coordinates, the condition reads 

The more elegant coordinate expression is obtained by lowering an index and reads 

■J^ = -2 {aik^ + 9jk^) ^ (62) 

which is indeed the defining relation for a special conformal Killing tensor as used in 
previous work (see e.g. jTHlEI)- An advantage of the present framework is that we do 
obtain an easy to handle and elegant, intrinsic expression also for the condition on the 
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tensor J in its type (1, 1) appearance, which is after all the way in which J is originally 
conceived. 

We finish this overview of the Riemannian case by briefiy re-deriving the most important 
properties of special conformal Killing tensors from ()54|1 . 

Theorem 2. If J is symmetric and satisfies \54\} for some function f G C°°(Q), then 
Nj = and / = tr J; moreover, if J is non-singular, then its cof actor tensor A is a 
Killing tensor. 

Proof. Acting with on (jSlI), and knowing that D^Xj = 0, it follows from the 
commutator property (j45|l that 

D5J= i(X®rf7-X/®D^0L). (63) 

We then get 

^"jxAY) - JmJ{Y)) - B"jyJ{X) + J(D?J(X)) 

The second term manifestly vanishes because d^di = d^d^L = 0. Taking X and Y to be 
basic for simplicity, the coefficient of Xf can be re-written as 

D-^(^^(X)) - B^jMY)) = D^jyigiT,X)) - D}'^(^?(T, F)) = giJY, X) - giJX,Y), 

which is zero in view of the symmetry of J. Lemma 2 now implies Nj = 0. 
From (jHH), we get 

V(trJ) = |((T,d«/)-(X;,^^z.)) 
= |(V/-(7(T,X^)) = V/. 

Hence, / = tr J (up to a constant, which is irrelevant). 

Finally, if A is the cofactor of J, meaning that J A = (det J)I, we have 

{D^J)A = -J{D'^A)+X"{detJ)I. (64) 

Again, it suffices to let X, Y, Z in what follows be basic vector fields, so that for example 
DxOl{Z) = D^{9l{Z)) = g{X, Z). For the sake of uniformity, we keep using the operators 
of the calculus along r, although everything here of course happens on the base space Q 
(and expressions such as X"{f) mean simply X(/)). From ()63|1 . it follows that 

D- J(Ar) = \ {{AYY{f)X - g{X, AY)Xf) . 

Using this to compute g{D'^J{AY), AZ) and taking a cyclic sum over X, Y, Z (indicated 
by an ordinary summation symbol), we readily obtain, knowing that also A is symmetric, 

J2 gmJiAY), AZ) = ^ (AXrif) g{Y, AZ). 
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Next, using this to compute ^ g{J{D^A)Y, AZ) = (det J) ^ g{D^A{Y), Z) via dSl, we 
arrive at 

(det J) J^gmMY), Z) = 5^(X-(det J) - ^^(r, AZ). 

But we know that Nj = imphes that dj{det J) = (det J) (i(tr J) (see e.g. E]), which 
can be written as rf^(det J)( JX') = (det J) rf^/(X'), for all X'. Taking X' = AX, it 
follows that 

X-(detJ) = (AXr(/), (65) 

which in turn leads to 

J29i^xAiY),Z) = 0. (66) 

This is the way to express that A, as type (1, 1) tensor, is a Killing tensor. The more 
familiar way is to look at the type (0, 2) tensor field A, obtained by lowering an index, so 
that JaA= (det J) g. The cyclic sum condition 

J]D«i(r,Z) = (67) 

then is equivalent to (jUBj) . □ 



6 An outlook for further study 

As stated in the introduction, our goal is to develop generalizations of the classical cases 
of Hamilton- Jacobi separable systems, or completely integrable systems, of which we have 
not given any examples here, because such examples can abundantly be found in the cited 
literature. There are reasons to believe that a tangent bundle approach will then have 
advantages over a cotangent bundle framework. The present study is a preliminary in- 
vestigation about understanding how things work on a tangent bundle. Even so, we have 
already obtained in Sections 3 and 4 some general results relating to an arbitrary La- 
grangian function (not necessarily one of 'mechanical type'). But for a full generalization, 
also more general type (1, 1) tensors J should be allowed, with components depending 
on coordinates and velocities (or coordinates and momenta). For such a J, the notion 
of complete lift is lost, so how to proceed? The point now is that we indeed have an 
idea of how to proceed in the tangent bundle set-up. It suffices to look at the expression 
(^Hj) for J'^ and to observe that the right-hand side is perfectly defined also for a tensor 
field J along the projection r. In fact, we are then talking about a more general lifting 
procedure, which has been fully developed already in fn\ |22] and is sometimes referred 
to as the F-lift. For a given J along r which is not basic, the formula 

JrJ = J" + (VJ)'', (68) 

indeed defines a type (1, 1) tensor field on TQ, which depends on a given second-order 
equation field F. A number of the calculations which will be involved in such a general- 
ization start off in exactly the same way as in the present paper, but of course without 
the simplifications coming from certain objects being basic. 
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A particular case of interest which could already significantly generalize the well-known 
Riemannian situation of the previous section, is to let g be the metric along r coming 
from a Lagrangian which is the square of a Finsler function F. It is then appropriate to 
start from a J along r (with the zero section of TQ excluded) which is homogeneous of 
degree zero in the velocities. Thus, one can use the lifting procedure above, where F is 
the Euler-Lagrange field of F^. Work along these lines is in progress. 

Appendix A: Aspects of integrability analysis 

Starting from an arbitrary J on Q and defining R on TQ by Q, the conditions for having 
CrR = are, according to Proposition 5, that J is symmetric and parallel, and further 
commutes with the Jacobi endomorphism $. But we have argued indirectly in Section 5 
that in the Riemannian case, the third condition must be an automatic consequence of 
the first two. An explicit verification of this fact can only come from an integrability 
analysis on the partial differential equations satisfied by J. 

The assumption is that VJ = 0, and since J is basic, also J = for all X, so that ()45j] 
implies that also D^J will be zero for all X (hardly a surprise of course in view of 
The next interesting commutator to look at here is [V,D^], or more generally [DJ,Dy]. 
Its general expression reads (see |22j ) 

Dy] = Bfx^Y]^ + D^(x,y) + /^Rie(x,y)- (69) 

Clearly, only the last term matters here; it is a derivation which, when acting on a vector 
field Z along r is given by 

fiRieix,Y)Z = Rie{XX)Z = -D^7^(X,r). 

In fact Rie is here simply the classical Riemann tensor. Recall also that we have the 
following relations linking $ and TZ, 

= 3n, $(X) = 7^(T, X), (70) 

which implies for example that $(T) = 0. 

Now, to express that [V,DJ]J must be zero, we compute 

[V,D5](jr) - J([V,D5]F) = Rie(T,X)(JF) - J(Rie(T, X)r). 

Using the second of (fTOj) and the property DyT = Y, we can write Dy7?,(T,X) = 
Dy$(X) + 71{X, Y), by which the above expression in the end reduces to 

A{X, Y) := -B'^jY^X) - 7^(X, JY) + J(D^$(X)) + J(7^(X, Y)) = 0. (71) 

In particular, knowing that $ here is quadratic in the velocities so that D,^$ = 2$, it 
follows that 

A{X, T) := -D;;^'^!^) - T^iX, JT) + J($(X)) = 0, (72) 
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from which we further obtain that 

= g{A{X, T), T) = -g{n{X, JT), T) - ^?(D:;t*(X), T) + J$(X), T). (73) 

Now we recall from the study of the inverse problem of the calculus of variations that also 
$ is symmetric with respect to g (see e.g. Theorem 8.1 in [22] )• It then follows from the 
Bianchi identity 

5^(7(7^(X,y),Z) = 0, (74) 
applied with arguments X, JY, T, that 

g{n{X, JY), T) = (?($ J(r), X) - (7(<f X, JY) = 0. (75) 

So the first term on the right in (j7H|) vanishes. Moreover, also is symmetric 

with respect to g for all X, which implies that the second term can be re- written as 
-c/(X,D5t'^'(T)) = ^(X,$(D>tT)) = c/(X,$J(T)). The conclusion from ^ therefore 
is that g{J^{X), T) = 0. Computing the Dy-derivative of this result and subtracting the 
same expression with X and Y interchanged, it follows by using the first property in (f7n|) 
and the symmetry of J and $ that 

3g{'JZiX,Y), JT) = (7((J$ - $J)X,F). (76) 

Next, we use the Bianchi identity again to write 

g{niX, Y), JT) = g{n{X, JT), Y) - g{n{Y, JT), X), 

and make use of (|72p to compute the right-hand side. Taking X and Y to be basic for 
simplicity, we can write (y'(Djrp$(X), Y) = Djrj,((7($X, Y)) and likewise for the term with 
X and Y interchanged. It then readily follows that g{TZ{X, F), JT) = 5f((J$-$J)X, Y). 
Comparison with ()76p leads to the conclusion that both sides must be zero, for arbitrary 
X, Y . Hence, we have shown in a direct way that in the Riemannian case, 

J = J and V J = ^ $ J = J$. (77) 

We had an indirect proof of this fact in Section 5. It is of interest to illustrate that it is 
indeed a non-trivial property by looking at coordinate expressions. 

We have 

$}=7^1^M^ 7^t^ = and thus = i^^^.wV, (78) 

where R\jj^ are the components of the Riemann tensor, and are skew-symmetric in the last 
two subscripts. Now, from Jj|^ = 0, taking a further covariant derivative, swapping indices 

and using the Ricci identities, the property which immediately follows is Jj-R^^; = R]rni'^i- 

But the commutation of J and $ is a different property and says that J]{Rkmi + ^Imk) ~ 
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Appendix B: Darboux-Nijenhuis coordinates 



It is well known that on a general (regular) Poisson-Nijenhuis manifold of dimension 
2n, if the recursion operator R has n distinct eigenvalues, their exist so-called Darboux- 
Nijenhuis coordinates, which diagonalize R and are at the same time Darboux coordinates 
for the symplectic form (see e.g. [SHI El)- This will apply in particular to the general 
situation on TQ, described in Sections 3 and 4. We wish to investigate here in some 
detail what the structure is of the eigenspaces of our R and how the construction of 
Darboux-Nijenhuis coordinates works when the eigenvalues are maximally distinct. 

We begin by establishing results which are valid without special assumptions on the type 
(1, 1) tensor J on Q, except that we will only consider real eigenvalues. 

Lemma 3. If ^ = X" + is an eigenvector of R, corresponding to the eigenvalue \, 
then 

JX = XX and UX + JY = XY. (79) 
It follows in particular that X is an eigenvector of J. 

Proof. Using the characterization of R as described by 1)3011311) . it is immediate to see 
that R^ = X^ is equivalent to the two relations ()79|) . □ 

Lemma 4. J and J have the same eigenvalues. In fact, if X is an eigenvector of J, then 
Xj g is an eigenform of J with the same eigenvalue. 

Proof. We have = g'-^Jlgij, and therefore 

l\-X5]=g'\Jl-X5l)g,,. 
Both statements now easily follow. □ 

Lemma 5. Suppose that J is non- degenerate and has n distinct eigenvalues (which then 
are non-zero). Then, if JX = XX , there exists a vector field Y along t, such that 
JY = XY - UX. 

Proof. From g{JX,Y) = g{X,JY) = Xg{X,Y), it follows that g{X,JY-XY) = 0,VF. 
Extending X = Xi to an orthogonal frame {Xi, . . . , X„} = {Xi, Xa} for g, and putting 
Y = a'^Xi, it follows that a^JXi — XXi) G sp {Xq}, Va*, which implies that 

JXi = XXi + btX^, z = l,...,n 

for some functions 6". We know that A is an eigenvalue of J as well, and that its eigenvalues 
are distinct. Hence, there exists a unique vector field of the form Xi + c"Xq, which spans 
the kernel of 7 - XI. But (7 - A/)(Xi + c"X„) = (6^ + c^6;^)X„, so the fact that 
unique functions exist which make this zero implies that det{b'^) ^ 0. Now consider 
the equation JY = XY — UX for the unknown Y = a^Xi G '^(t). Since g{UX,Y) is 
skew-symmetric in X,Y, we know that g{UXi, Xi) = and thus UXi = d'^Xa for some 
functions c?". The equation for Y can now be written in the form a^b'^ = — — a^6" and 
clearly has a unique solution for the for each arbitrary choice of a^. □ 
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Proposition 9. Let J be diagonalizahle with distinct non-zero eigenvalues. Then a com- 
plete set of eigenvectors of R can he constructed as follows: (i) let Xi denote the eigen- 
vector of J with eigenvalue Aj and Zi the eigenvector of J with the same eigenvalue; (ii) 
for each Xi, construct a vector Yi such that JYi = XiYi — UXi. Then Zi^ and Xi" + Yi^ 
are eigenvectors of R, corresponding to the eigenvalue Aj. 

Proof. We have 

R{Zi^) = {JZi) = XiZi^ , 

Rix," + Yn = (jx,)" + (f/x,)" + (7ri)" = A,(x,- + Fn, 

from which the result follows. □ 

This much for the purely algebraic aspects. Now let us further assume that Nj = 0. 
Darboux-Nijenhuis coordinates in fact should do three things at the same time: not only 
diagonalize R in coordinates, but also separate it, and bring the symplectic form into 
canonical form. It was proved in the fundamental paper of Frolicher and Nijenhuis JH] 
that if J is (algebraically) diagonalizable and the eigenvalues have constant multiplicity, 
then the necessary and sufficient condition for diagonalizability in coordinates is that 
TCj = 0, where the Haantjes tensor TCj can be defined by 

HjiX, Y) = J'Nj{X, Y) + Nj{JX, JY) - JNj{JX, Y) - JNj{X, JY). (80) 

Obviously, Nj = implies Hj = 0, but evaluating Nj on eigenvectors X and Y belonging 
to different eigenvalues. A, /i say, further gives 

= NjiX, F) = (A - fi){X{fi)Y + r(A)X), 

so that X{fi) = Y{X) = 0. Hence, in coordinates which diagonalize J, the eigenvalues 
will only depend on the coordinates of the corresponding eigendistribution, which is the 
meaning of saying that J is separable in coordinates. Conversely, if J is separable, one 
can verify in such coordinates that Nj = 0. In other words, Nj = (for a J which has the 
algebraic properties stated above) is the necessary and sufficient condition for separability 
in coordinates. Note in passing that the tools for studying such issues when J would more 
generally be a tensor field along r have been developed in |23j . 

To understand what happens with R on TQ now, we need to look at the expression of R 
in a coordinate basis, rather than in the adapted frame as in ([HUjl : it reads (still for the 
general situation described by Proposition 1) 

R = J]-^ ® dq^ du^ + iU] + - J-Ti)^ (g) dq^. (81) 

The following procedure now will lead to Darboux-Nijenhuis coordinates. First perform 
the Legendre transform (g, m) {<1,P = dL/du). Even though this is to be regarded 
here as a change of coordinates on TQ, the result will be that R acquires the form of the 
complete lift J on T*Q, i.e. 

R = J] [-Q-®dq^ + Q^^^dpj +p, I ^ - ^ 1 —®dq^. (82) 
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The 2-form ujl meanwhile will already take its canonical form in the (a;,p) coordinates. 
Now, assuming that J has distinct eigenvalues and zero Nijenhuis torsion, we know that 
there exists a coordinate change on Q which will diagonalize the expression Jj{d/dq'^®dq^) 
in such a way that the eigenvalues depend on at most one new coordinate q'. The resulting 
point transformation (g, u) — > (g', u') on TQ, when expressed in the non-tangent bundle 
variables {q,p), formally is a 'canonical transformation' {q,p) — >■ {q',p'), i.e. it defines 
another Darboux chart for the symplectic form ljl and it will have the additional effect 
of diagonalizing (and separating) R. 

From a tangent bundle point of view, the first step in this procedure is rather unnatural, 
because it is not a tangent bundle change of coordinates. At first sight, it may look like 
one should nevertheless not change the order of the operations, because even though J 
and J have the same eigenvalues, a coordinate transformation which diagonalizes J will 
generally not at the same time diagonalize J. However, the two coordinate changes under 
consideration here are of course of a quite special type: a Legendre transformation which 
does not change q but changes the fibre coordinate, and a point transformation. It is 
clear that such coordinate changes commute, so one can just as well diagonalize J first 
and then the subsequent Legendre transform will not destroy the diagonal form of J, will 
bring cul in canonical form, and at the same time will take care of the diagonalization of 
J. 

That the reversed procedure is somewhat more natural for the tangent bundle set-up may 
become clear in the special case that J is symmetric. It then follows from gijJl = QikJ] 
that in coordinates which diagonalize J, we will have gkj{\'^^) — A^-^-*) = and thus 
Qkj = for j ^ k. This gives useful information also when there is no urge to pass 
to Darboux-Nijenhuis coordinates: it means that in coordinates which diagonalize J, the 
given Lagrangian will separate with respect to the velocity variables, i.e. it will become 
of the form: L = L\q, m*), where L* depends on only. 
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